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Abstract 

A pairing model for nucleons, introduced by Richardson in 1966, which describes 
proton-neutron pairing as well as proton-proton and neutron-neutron pairing, is re- 
examined in the context of the Quantum Inverse Scattering Method. Specifically, 
this shows that the model is integrable by enabling the explicit construction of the 
conserved operators. We determine the eigenvalues of these operators in terms of 
the Bethe ansatz, which in turn leads to an expression for the energy eigenvalues of 
the Hamiltonian. 



1. Introduction. 



Pairing model Hamiltonians have again become the focus of many theoretical con- 
densed matter investigations due to the fact that experimental work on metallic nanopar- 
ticles ( also refered to as small metallic grains ) has detected evidence of pairing interac- 
tions In order to gain an insight into the physical properties of small metallic grains, 
substantial attention has been devoted to the analysis of the ( reduced ) BCS Hamiltonian 
which is beheved to be appropriate to describe the dynamics of these systems . An im- 
portant point in this regard is that the treatment originally proposed by Bardeen, Cooper 
and Schrieffer for bulk systems, using variational wavefunctions with an undetermined 
number of particles ( grand canonical ensemble ), is not applicable to the study of a small 
metallic grain where the number of electrons remains fixed ( canonical ensemble ). This 
aspect has generated activity in analysing the BCS Hamiltonian under this constraint 
0. Remarkably, an exact solution of the reduced BCS Hamitonian was obtained some 
time ago in a series of works by Richardson and Sherman [Q. In these papers, the BCS 
Hamiltonian was studied for the purpose of application to pairing interactions in nuclear 
systems, and as such the solution escaped the attention of the condensed matter physics 
community for a considerable time. More recently, it was shown by Cambiaggio, Rivas 
and Saraceno that the model is also integrable; i.e. there exists a set of mutually com- 
muting operators which commute with the Hamitonian. These features can be reproduced 
in the framework of the Quantum Inverse Scattering Method ( QISM ) using a solution of 
the Yang-Baxter equation as shown in [^, 0] . This approach has the significant advantage 
that the computation of form factors and correlation functions can be undertaken in this 
algebraic framework R M. 



In the papers |y, |T0| Richardson introduced a coupled pairing model for nuclear systems 
which accomodates proton-neutron pairing interactions as well as the proton-proton and 
neutron-neutron couplings. This model was recently studied in ||TT|. In this paper we 
will show that the techniques employed in H, 0] can be applied to this model to show 
that it is integrable, and for the determination of the energy spectrum. This formulation 
also opens the possibility for the calculation of form factors and correlation functions by 
algebraic means, in analogy with the results of ||^. For the present case the solution is 
obtained through the use of a solution of the Yang-Baxter equation associated with the 
Lie algebra so(5). The Hamiltonian has the explicit form 

D D 

H = Y. '^""^ - 9 E (^Kl)^'^(l) + &3(2)&fc(2) + &](3)&.(3)) (1) 

where g is an arbitrary coupling parameter and D is the total number of distinct en- 
ergy levels. Above, rij is the number operator for paired nucleons at energy level ej 
and b{i)j, b^j{i), i = 1,2,3 are the annihilation and creation operators for three sets of 
generalized, non-commuting, hard core boson operators satisfying the relations (amongst 
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others) 



{b]{z)y = 0, t = l,2, (&](3))' = 
[b,{^, hil)] = [b]{i), blil)] = [6,(0, blil)] = 0, for k ^ J. 



The three sets of hard core boson operators correspond to proton-proton, neutron-neutron 
and proton-neutron pairing in a nuclear system. Their exphcit forms will now be made 
clear. Two sets of two-fold degenerate Fermi operators c± and d± and their hermitian 
conjugates are introduced. c± and d± represent the protons and neutrons, the subscripts 
± referring to time reversed states. The hard-core bosons are realized by 



with appropriate definitions for the hermitian conjugates. The Hamiltonian (|l]) corre- 
sponds to a special case where the energy level spacings for the protons and neutrons are 
the same and the scattering coupling is the same for proton-proton, neutron-neutron and 
proton-neutron pairings. In this instance the Hamiltonian acquires an isospin symmetry 
which plays an important role in our analysis below. 

The energy levels are degenerate. Each level can be empty or occupied by protons 
and/or neutrons in two-fold degenerate time reversed states. This gives the degeneracy of 
each level as 2^ = 16. However, the above Hamiltonian only scatters paired nucleons giving 
rise to a blocking effect ( c.f., the blocking effect for the BCS model discussed in |0 ). For 
any energy level which contains an odd number of nucleons, the pairing interaction acts 
trivially, and these states can be discarded from the Hilbert space. ( Hence, we choose 
Uj to count the number of paired nucleons at e,, rather than the number of nucleons. ) 
Furthermore, it is assumed that the proton-neutron pairing is between time reversed states 
which are symmetric under interchange of protons and neutrons. In this case the pairing 
interaction is also trivial on the non-time reversed paired proton-neutron states and the 
antisymmetric time reversed paired proton-neutron state, giving rise to a 5-dimensional 
space at each level on which the scattering is non-trivial. ( In the language of |0, 
this Hilbert space is spanned by the Seniority- Zero states. In the subsequent paper 
Richardson extended his results to Seniority-One and -Two states. However, we will not 
consider these cases here. ) It is convenient to use the fundamental representation of the 
so(5) Lie algebra to construct the local operators acting on each of these spaces, which 
we will now discuss. 



We can construct the fundamental representation of the so(5) Lie algebra in the fol- 



b\l) 
b\2) 

6t(3) 




{ cidi + did 




2. The Lie algebra so(5). 
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lowing manner. Define m = 6 — m. For 5x5 matrices, consider the subset 



'^m 



where denotes the matrix with 1 in the (m, n) position and zeroes elsewhere. Note 
that = 0. We will denote the 5-dimensional space on which these operators act by 
V . The operators close to form the fundamental or defining representation of the Lie 
algebra so(5) with commutation relations 

K, af] = Car + 4«; + 5joFm + (2) 
A basis for the Lie algebra is given by the set 

{a™ : 1 < m < n < 5} 
which gives ten linearly independent generators. Explicitly the basis generators read 



a\ 


= e\ 


-el 
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-el 
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4 
«5 


= e\ 


-e\. 



Note that the representation is unitary, and specifically 

(«n )^ = (3) 

Next, we recall some established results on the representation theory of so(5). For a 
more detailed discussion, see for example |T^. Finite dimensional irreducible representa- 



tions of so(5) are uniquely determined by the highest weight labels Ai, A2 which are the 
eigenvalues of the Cartan elements 

h\ = ctg, hi = CI4 

acting on the highest weight state. ( These operators are self-adjoint and mutually com- 
muting and so can be diagonalized simultaneously. ) The highest weight state is the 
unique vector v which vanishes under the action of the raising operators; viz, 

= a^^v = alv = alv = 0. 

The highest weight labels Ai, A2 take integer or half odd-integer values and are subject 
to the conditions 

Al > A2 > 0, Al - A2 G Z. 

In the case of the fundamental representation we have Ai = 1, A2 = and the highest 
weight vector is |1) = (1, 0, 0, 0, 0)*, with t the matrix transposition operation. 
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The so(5) algebra admits a second order Casimir invariant 

5 

C = Y,a^al (4) 

m,n 

commuting with all elements of so(5), which can be verified explicitly from the com- 
mutation relations (0). Due to Schurs lemma, on each irreducible finite dimensional 
representation the Casimir element (^) takes a constant eigenvalue, which is given by 

Xc(Ai, A2) = 2 (Ai(Ai + 3) + A2(A2 + 1)) • 

For the fundamental representation we clearly have Xc(l? 0) = 8. 

An important ingredient in the following construction is the existence of a canonical 
so(3) subalgebra spanned by the basis elements 

I-P — = L = a"^. 

The Casimir operator for this so(3) subalgebra is given by 

K = L+L- + L-L+ + {Ly. 

The irreducible finite dimensional representations of the so{3) algebra have a unique 
highest weight vector w which satisfies L~^w = 0. These representations are uniquely 
characterized by the eigenvalue fi of L° acting on w. The allowable values of /i are such 
that it is a non-negative integer or half odd-integer. The eigenvalue of K on such a 
representation is given by 

The realization of this canonical so(3) subalgebra we give below is referred to as the 
isospin algebra in 0. 



3. The Yang-Baxter equation and integrability. 



The basis for constructing an integrable model through the QISM ||T^ is a solution of 
the Yang-Baxter equation 

Ri2{u - v)Rn{u)R23{v) = R23{v)Ri3iu)Ri2{u - v) (5) 

which is a matrix solution acting on a three-fold tensor product space V ®V ®V . Above, 
the subscripts refer to which two of the three spaces the operator R{u) G End(y ® V) 
acts upon. Solutions of the Yang-Baxter equation associated with representations of 
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Lie algebras are well known. The i?-matrix solution associated with the fundamental 



representation of so(5) discussed above takes the following form |jT^: set 

5 
5 

m,n 



en ® e^. 



m,n 



Then 



Riu)=I + ^P-^Q 

U U + 67] 



provides a solution of (|]) with rj being a free parameter. This solution is so(5) invariant 
in that 

[R{u), / (g)x + x(g)/]=0 

for any x G so(5). We note the properties: 
1. Unitarity 



R{u)R{-u) = { ^^—^ 



2. Crossing symmetry 



- ^V) = {1® A)R{u){I ® A) 

where A is the matrix with elements A^ = 6^ and ti denotes matrix transposition in the 
first space of the tensor product. 

By the usual procedure of the QISM we define a transfer matrix acting on the D-fold 
tensor product space via 

t{u) = tro {GqRqd{u - to) Rm{u - ti)) 

which gives a commuting family satisfying \t{u), t(f )] = 0. Above, tr^ denotes the trace 
taken over an auxiliary space labelled by and G can be any matrix which satisfies 

[R{u), G] = 0. 

We choose G = explarjal) for this particular model. Using either the analytic Bethe 
ansatz, which exploits the unitarity and crossing symmetry properties [|15|, or the algebraic 
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Bethe ansatz developed by Martins and Ramos for the so{n) series [|rB|, the eigenvalues 
of the transfer matrix are found to be 



AM = exp(c,)n '"r'*^f"' n i" ""•"'! 

^ V (« - f-t ) V (« - V, + >)) 

+ -p(-"") n (.-,,+3,,) n („.-..+2„) + ''"("^ 

where Ao('u) are the transfer matrix eigenvalues for the i?-matrix associated with the spin 
1 Babujian-Tahktajan model [O], with inhomogeneities fj. These eigenvalues read 



(m - f j + 7]) ^ -^ {u - Wj + 2r]) 

I TT ju-V,) -TT - + 3^) 

J-.J- (m — t>j + 2rj) (li — + 7]) 
{u — Wj + 3?7) {u — Wj) 



N 



+n 



u — Wj + rj) {u — Wj + 2rf) 



The parameters Vj, are required to satisfy the Bethe ansatz equations 

expH)!!-^ r = -\\t^ ^tIIt^ 

Y (^J' ~ ~ ^> I ^'"i -'^i + V) i [Vj -Vi- 2r]) 

TT (wj -Vj-y) ^ TT jwj -Wk-y) 
i i^j -Vi + v) \^ {wj -Wk + v)' 



T, = hm (6) 

u^f-j 27] 



Define the operators 
which satisfy 

T,] = 0. (7) 

Explicitly 

Tj = GiRjoi^j - eD)---Rj{j+i){ej - ej+i)i?j(j_i)(ej - ej_i)...Rji{ej - ei). 
Now, by taking the quasi- classical expansion 

T^ = I + ^r^+oiv') (8) 

we find 

(pjk 



Tj = a 



where 

5 

m,n 

and for ease of notation we set ■?/' = Og. Setting 

5 

m,n 
5 

5 

m,n m,n 
5 

m,n 

= 20 

shows that we may write 




which satisfy [tj, t^] = in view of 

Consider the action of the so(3) Casimir on the D-fold tensor product space 

K^f^iLtLj + LrLl + U^L^). 
hi 

It is easily deduced that 

[K, = 0. 

When a = { and so G = / ), 

[K, T,] = 

since in this instance the operators tj are so(5) invariant as a consequence of the so(5) 
invariance of the i?-matrix. We thus see that 

[K, r,] = 

in general. The set of operators {rj, K} are mutually commutative and so can be used 
to define an integrable Hamiltonian through any function of these operators. With an 
appropriate choice, we show below that the pairing Hamiltonian (0) introduced in |^ can 
be reproduced, thus establishing the integrability of this model. 
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4. Pairing Hamiltonian. 



First, let us realize the 5-dimensional space V in terms of the two sets of two-fold 
degenerate Fermi operators c±, Cj_ and d±, dj. introduced earlier. We make the identifi- 
cations 



1) 
2) 

3) 

4) 
5) 



|0), 

dldl |0) , 

-^[cldl + dlcl) |0) 

del |0) , 
clcUldl |0) . 



Set n'^ = c^^c+ + cLc_, n'^ = dl^d^ + d^_d_ and n = l/2{rf + n''), which measures the 
number of paired fermions. We have the following realisation of the so(5) generators 

t J 



c_c 



dldl, 
72 V 



72 



1 ( c 

- (n 

2 ^ 



2 ^ 



and 



.2U 



.3U 



The representation of the canonical so(3) subalgebra generated by {L^, , L^} is the 
isospin algebra of and the operator ip is a. U{1) generator. We can also identify the 
generalized hard core boson operators with certain elements of the so{5) algebra through 



6t(l) 



1, b\2)=al b\3) 



and corresponding relations for the hermitian conjugates from 
1 



We may now express 



':e = 6^(1) ® 6(1) + +61"(2)® 6(2) +6t(3) 6(3) 

+6(1) (g) 6^(1) + 6(2) (g) b^{2) + 6(3) ® 6^(3) 

+L+ (g) + (g) + L'^ (g) + ijj (g) ^p. 
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Define the Hamiltonian 



H 



^ / D g \ 1 ^ ^ 

7 E ( - - ) - ^ + 4 D + ^ E -.-^ + E 

\ 3 / 0,k j 

D D D \ ( ^ \ ^ 

E(^.- - -)^^- - ^ E^^-^ + - E^^-^^ + ^ + E^^- 



where in the last line we have use the fact that the so(5) Casimir invariant C takes the 
eigenvalue 8 in the 5-dimensional representation, as mentioned earlier. Expressing the 
so(5) elements in terms of the hard core boson operators as indicated above yields (|l|) 
with g = 2/ a. The Hamiltonian describes two coupled identical BCS models, where, in 
addition to the customary pairing ( characterized by the operators &^(1), &^(2), 6(2) ), 
fermions from each BCS system at the same energy level and in time reversed states can 
pair in such a fashion that the wave function is symmetric under interchange of the two 
BCS systems ( described by 6^(3), 6(3) ). This shows that the Hamiltonian has a natural 
interpretation as a pairing model for nucleons which includes proton-neutron pairing. 



5. Energy Spectrum 



It remains to determine the eigenvalues of the Hamiltonian, which is achieved by 
computing the eigenvalues of the conserved operators. The eigenstates can be labelled by 
the eigenvalues of the Cartan elements hi = a\ = h2 = a\ = L° acting on the tensor 
product space. For each value of M and which appear in the solution of the Bethe 
ansatz equations, we find that the corresponding eigenvalues of the Cartan elements are 

Ai = D-M, A2 = M-N = ^ 

or equivalently 

n = M, n'^ -n^ = 2{M - N). 

We begin with the operator K. It follows from the algebraic construction of the Bethe 
states due to Martins and Ramos |16| that each Bethe state is a highest weight state with 
respect to the so(3) subalgebra. ( When a = the Bethe states are highest weight states 
with respect to the full so(5) algebra, but generic values of a break this symmetry. ) 
Since K is simply the so(3) Casimir operator, it takes the eigenvalue Xk{M ~ N) = 
{M - N){M - N + 1) on such a Bethe state. 

From (P,||) we see that that the eigenvalues for Tj can be obtained from the quasi- 
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classical limit of the transfer matrix eigenvalues. These read 

o 2 ^2 
A ■ = a + V V 

and the Bethe ansatz equations take the form 

^2 ^4 ^2 
«+> = > +> , 

V ~ 1^3 / ^' ~ 

M , N 



E— = E— 



Using these Bethe ansatz equations we can derive the following identities 

N M ^ N N 

EE^— = EE^^ 

M D M M . M N 

aM+vv^- = EE^+EE^— 

= 0, 

M N M N 

3 I ^ ^ 3 I ^ ^ 

N M N N 



M D M D 

M M D „ MM. M N ^ 

«E-. + EE^^ - EE^+EE^^ 

J 3 k-' J ijLj -I 3 I 

= 2M{M -1)-2MN + N{N -1). 

We now obtain from these identities 

D 

Y,\, = a{D-M), 
j 

D D M 

^ejXj = a^ej + D{D-l)-a^Vj + 2M{M -l) + N{N -1)-2M{D + 

j 3 j 
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Using these results we find that the energies are given by 



E 



a 

M 



/ D \ 1 ^ ^ 

\ j / j,k j 



This energy expression shows that the roots {vj} of the Bethe ansatz equations are simply 
the quasi-particle excitation energies. It is interesting to note in the case = ^ n'^ = 0, 
where the model describes a single reduced BCS system since there is only one type of 
nucleon, the Bethe ansatz equations and energy expression coincide with those obtained 
by Richardson and Sherman 0]. 

Finally, it is necessary to compare our results with those obtained by Richardson |P, 
in which the following Bethe ansatz equations were obtained 

1 A 1 _ M(M-3) + (M- A^)(M-iV + 1) y> 1 

and are obviously different from our results. The explanation for this difference stems 
from the fact that the ansatze adopted for the eigenstate wavefunctions are different in 
each case. Richardson chose wavefunctions which have eigenvalue zero under the action of 
the isospin operator LP. As we have indicated above, the states we construct are highest 
weight states with respect to the isospin algebra. An important open problem is to prove 
the equivalence of these two solutions. 



6. Conclusion. 



We have shown, by using the QISM, that the coupled BCS Hamiltonian proposed by 
Richardson |^ to accomodate proton-neutron pairing in nuclear systems is integrable. We 
have also determined expressions for the energy eigenvalues of the model in terms of a 
Bethe ansatz solution of coupled equations. It should be emphasized that although the 
model studied here is based on a specific Lie algebra and representation, the construction 
that we have employed to demonstrate integrability is entirely general. It can be equally 
applied to any representation of any Lie algebra or superalgebra, to yield a vast class of 
integrable systems. For a recent example based on the spin 1 representation of the so(3) 



algebra see |21 



An interesting question to ask is whether this solution is complete; i.e., are all energy 
levels obtained? It is well known that for many Bethe ansatz solvable models where 
there is an underlying Lie algebra symmetry the eigenstates are highest weight states 
with respect to this algebra O, EO]. By computing the dimensions of each multiplet 



12 



generated by these highest weight states and then employing a combinatorial argument, 
completeness can be proved. 

For the present model, where the i?-matrix solution has so(5) symmetry, this symmetry 
is broken in the construction of the transfer matrix by the inclusion of the operator G, and 
an so(3) symmetry ( isospin ) remains for the conserved operators Tj. The degeneracies of 
the eigenvalues can be classified in terms of so(3) multiplets. In the a — >■ limit the so(5) 
symmetry is restored and we find that the Bethe states are so{5) highest weight states, 
so there is an increase in the degeneracy of each eigenvalue at this point, or equivalently, 
a decrease in the number of distinct eigenvalues. Fortunately, the Bethe ansatz equations 
(|lO) admit more solutions for generic values of a than the a = case and automatically 
accomodate this facet. This is easily illustrated in the instance D = 2,M = 1,A^ = 0, in 
which case we need only solve 

a + ^— + = 0. (12) 

V - ei V - €2 

For non-zero a, this is a quadratic equation with two solutions for v. When a is zero, the 
equation is linear with the unique finite solution f = (ei + e2)/2. ( The equation is also 
satisfied hj v = oo which is the limit of one solution of (|12D as a ^ 0. However, such 
infinite solutions are trivial in the sense that they do not contribute to the eigenvalues 
(§) for the conserved operators. ) For general values of D with M = 1, = 0, (^^) gives 
rise to a polynomial equation of order D ioi a ^ 0, but this equation reduces to order 
(D — 1) when a = 0. Whether the Bethe ansatz solutions give the complete spectrum for 
generic a is an open problem still to be solved, but the discussion above shows that it is 
possible since the breaking of the so(5) symmetry to so(3) is accompanied by an increase 
in the number of solutions to the Bethe ansatz equations. For the purpose of counting 



the states, the results discussed in |^ may be appropriate. 

A final aspect for consideration is the possibility to compute form factors and correla- 
tion functions for this model. By rederiving the solution in the framework of the QISM, 
we hope to motivate further studies that are necessary for this task, such as an analogue of 
Slavnov's formula for wavefunction scalar products, which is well known for su{2) models 
2|. 
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